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ABSTRACT 

Until recently, many of the dozens of quantitative predictions of the ambipolar-diffusion theory of 
gravitational fragmentation (or core formation) of molecular clouds have been confirmed by obser- 
vations and, just as importantly, no prediction has been contradicted by any observation. A recent 
paper, however, claims that measurements of the variation of the mass-to-flux ratio from envelopes to 
cores in four clouds decreases, in direct contrast to a prediction of the theory but in agreement with 
turbulent fragmentation (in the absence of gravity) and, therefore, the ambipolar-diffusion theory is 
invalid (Crutcher et al. 2008). The paper treats magnetic-field nondetections as if they were detec- 
tions. We show that the analysis of the data is fundamentally flawed and, moreover, the comparison 
with the theoretical prediction ignores major geometrical effects, suggested by the data themselves 
if taken at face value. The magnetic fluxes of the envelopes are also miscalculated. We carry out a 
proper error analysis and treatment of the nondetections and we show that the claimed measurement 
of the variation of the mass-to-flux ratio from envelopes to cores is not valid, no contradiction with 
the ambipolar-diffusion theory can be concluded, and no theory can be tested on the basis of these 
data. 

Subject headings: ambipolar diffusion - ISM: clouds, evolution, magnetic fields - MHD - polarization 
- stars: formation - turbulence 



I. INTRODUCTION 

The ratio of the mass and magnetic flux of interstel- 
lar molecular clouds has received well-deserved obser- 
vational attention in recent years (e.g., Crutcher 1999; 
Heiles & Cructher 2005). For a cloud as a whole, the 
mass-to-flux ratio is important input to the ambipolar- 
diffusion theory of fragmentation (or core formation) 
in molecular clouds. For fragments (or cores) within 
a cloud, the ambipolar-diffusion theory makes specific 
predictions of the mass-to-flux ratio and its evolution 
in time. Until now, observations have been in excellent 
quantitative agreement with the theoretical predictions 
(e.g., see Crutcher et al. 1994; Crutcher 1999 and correc- 
tion by Shu et al. 1999, pp. 196 - 198; Ciolek & Basu 
2000). Recently, however, Crutcher et al. (2008) (here- 
inafter CHT) have looked at the variation of the mass-to- 
flux ratio from the envelopes to the cores of four molecu- 
lar clouds and concluded that there is a fundamental dis- 
crepancy between the observations and the predictions of 
the ambipolar-diffusion theory. 

In what follows, we examine critically the claimed dis- 
crepancy. Before we undertake a direct re-examination 
of the data and the conclusions drawn from them, it is 
important that we summarize the key elements of the 
ambipolar-diffusion theory of protostar formation. This 
is especially necessary because, often, what is a predic- 
tion of the theory is quoted as if it were a requirement 
or assumption, and vice versa, i.e., input of the theory 
taken directly from observations is often quoted as if it 
were a prediction. The ambipolar-diffusion timescale is 
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also often a frequent victim of abuse in the astronomical 
literature. In this context, it is of utmost importance to 
clarify the difference in the way the process of ambipo- 
lar diffusion affects the ordered (or mean) component of 
the magnetic field in a molecular cloud and the magnc- 
tohydrodynamic (MHD) waves (or turbulence) that are 
superimposed on the mean field. It is also useful, for 
putting the current controversy in its proper context and 
perspective, to summarize how theory and observations 
of interstellar magnetic fields have interacted over the 
years and what the outcome of that interaction has been. 

We summarize in ii2.1l the main elements of the 
ambipolar-diffusion theory in a context relevant to the 
present controversy and to maintaining proper perspec- 
tive. In §2. 21 we first summarize the observational proce- 
dure followed by CHT and state their main conclusion, 
and then we state three major flaws in their analysis as 
well as three additional such flaws in their comparison of 
the results with the ambipolar-diffusion theory. 

In fj3l we explain the proper procedure for combining 
multiple observations of the envelope magnetic field in 
each cloud to derive an estimate of the mean and an 
associated uncertainty. An explicit assumption of the 
analysis by CHT is that line-of-sight effects (as encoded 
in cos6') are the same for the core and the envelope. An 
implicit assumption of the same analysis is that the mag- 
netic field does not vary intrinsically from one envelope 
location to another. None of these assumptions is jus- 
tified by the observations. Relaxing them in a manner 
suggested by the data themselves alters the manner in 
which errors should be treated and increases significantly 
the uncertainties of all results. 

In m we discuss the difference between a detection and 
an upper limit of the mean envelope magnetic field. Most 
of the individual observations yielded nondetections, a 
fact recognized by the authors. The authors state that 



for this experiment detections are not necessary. No up- 
per limits are calculated, which would be appropriate 
for nondetections; instead, nondetections and detections 
are treated in the same way with respect to error prop- 
agation and handling, as well as for quoting results. We 
remedy this problem and explicitly calculate upper limits 
where appropriate, using a full Monte Carlo analysis of 
the error propagation to properly account for the (large) 
errors propagating through a nonlinear equation (see eq. 
[7] below). 

In ^ we elaborate certain theoretical issues (points 
3a, b, c of i 32.2p which should be paid due consideration 
if a meaningful comparison between theory and observa- 
tions is to be achieved at least as far as the variation of 
the mass-to-flux ratio from the envelope to the core of a 
molecular cloud is concerned. 

2. THE DEVELOPMENT OF THE AMBIPOLAR-DIFFUSION 
THEORY. ITS MAIN ELEMENTS, AND RECENT OH 
ZEEMAN OBSERVATIONS 

2.1. Key Elements of the Ambipolar- Diffusion Theory 
of Protostar Formation 

Ambipolar diffusion (the motion of the two poles of 
charge, ions and electrons and possibly charged grains, 
which at densities characteristic of molecular clouds are 
well attached to magnetic field lines, relative to the neu- 
tral particles) is an unavoidable process in weakly ion- 
ized, self-gravitating systems, e.g., molecular clouds. It 
occurs on a timescale given by (see Mouschovias 1979b) 



TAD = 1.8 X 10*^ 



(tf^) 



yr, 



(1) 



where Xi is the degree of ionization. Unlike the early un- 
derstanding of this process, which envisioned magnetic 
flux escaping a cloud as a whole and thereby depriving a 
cloud of magnetic support and allowing it to collapse and 
form stars (e.g., Mestel & Spitzer 1956; Nakano & Tade- 
maru 1972; Nakano 1973, 1976, 1977), a re-examination 
of the process concluded that ambipolar diffusion is much 
more effective in the deep interiors of molecular clouds, 
where the degree of ionization is smaller than 10^^, than 
in their envelopes, where the degree of ionization is usu- 
ally greater or much greater than 10~^. The impor- 
tant difference in the two pictures is that, in the lat- 
ter, a cloud as a whole does not lose any magnetic flux. 
Instead, neutral particles, under their own self-gravity, 
contract slightly more rapidly than the plasma (and the 
field lines), thereby redistributing the mass in the cen- 
tral parts of a dense cloud and increasing the mass-to-flux 
ratio in the central flux tubes (but not necessarily in an 
individual fragment or core; see below). 

Early efforts to understand the possible role of am- 
bipolar diffusion in cloud dynamics and star formation 
consisted almost exclusively of comparing tad (actually, 
an earlier version of this timescale obtained by Spitzer 
1968 for uniform, cylindrical cloud models and which was 
greater than that given by cq. [1] by almost a factor of 3) 
and the free-fall time, which for a spherical cloud model 
of uniform neutral density p^ is given by 

1/2 /in3„_-3\l/2 
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Since tad was invariably found to be greater than Tff, it 
was always concluded that ambipolar diffusion is not a 



relevant process in star formation. The presumption was 
that a cloud gives birth to stars on the free-fall timescale, 
as orginally suggested by Hoyle in his hierarchical frag- 
mentation picture, or on a somewhat longer (" 'magneti- 
cally diluted" ) free-fall time because of the presence of a 
magnetic field frozen in the matter (Mestel 1965). 

A series of papers initiated and developed the main el- 
ements of the modern theory of single-stage fragmenta- 
tion (or, core formation) and star formation in molecular 
clouds as a natural consequence of primarily magnetic 
support of such clouds (Mouschovias 1977, 1978, 1979, 
1981, 1987a,b, 1989, 1991a, b, c; Shu et al. 1987). Obser- 
vations were already suggesting that star formation does 
not occur on the free-fall timescale but on a timescale 
exceeding 10^ yr (e.g., see Roberts 1969; Zuckerman & 
Palmer 1974). Magnetic support was neither an arbi- 
trary theoretical assumption nor an ingenious theoretical 
conclusion. It was a simple consequence of the fact that 
the then observed interstellar, large-scale magnetic field 
of 3 pG can easily support a mass of 5 x 10^ M© at the 
approximate mean density of the interstellar medium of 
lcm~^. Moreover, atomic-hydrogen clouds, thought to 
be the progenitors of molecular clouds have never been 
observed to have mass-to-flux ratios exceeding the criti- 
cal value for collapse 
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obtained from nonlinear, exact equilibrium calculations 
(Mouschovias & Spitzer 1976). (The quantity G is the 
universal gravitational constant.) In addition, if clouds 
were to begin their lifetimes in a magnetically supercrit- 
ical state, contraction (ordered) velocities characteristic 
of collapse (several kms~^) would be commonplace. In- 
stead, the gravitational potential energy of typical HI 
clouds is smaller than their magnetic energy by two or- 
ders of magnitude, and velocities characteristic of col- 
lapse for a cloud as a whole are not observed even for 
molecular clouds. 

Mouschovias (1987a) described a new theory of star 
formation in primarily magnetically supported molecu- 
lar clouds, with hydromagnetic waves (or turbulence) 
contributing part of the support against gravity. (At 
the time, it was referred to as a "scenario" because the 
rigorous, quantitative calculations needed to justify the 
term "theory", which implies predictive power, had not 
yet been undertaken.) Ambipolar diffusion, which as ex- 
plained above is an unavoidable process in a weakly ion- 
ized, self-gravitating physical system such as a molecular 
cloud, acts differently on the large-scale (or mean) mag- 
netic field and on short-wavelength disturbances associ- 
ated with MHD waves or turbulence. Because of am- 
bipolar diffusion, disturbances of wavelength A smaller 
than the Alfven lengthscale 



A a = TTWATn 



(4) 



cannot be sustained in the neutrals; i.e., they decay 
because of neutral-ion friction. (The quantity vx — 
i3/(47r/9n)^/^ is the Alfven speed in the neutrals, and T„i 
the mean (elastic) collision time of a neutral particle in 
a sea of ions.) Unlike the timescale for redistribution of 
mass in the central flux tube of a cloud given by equation 
II]), the timescale for the decay of such disturbances is 



wavelength (A) dependent and is given by 
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. (5) 

This is a very short timescale even relative to free-fall, 
and it becomes smaller the smaller the wavelengths one 
considers. For typical molecular-cloud parameters, the 
Alfven lengthscale is arithmetically equal (0.3 pc) to 
the thermal critical lengthscale (essentially the Bonnor- 
Ebert radius for collapse of an isothermal sphere sup- 
ported by thermal pressure against its self-gravity). 



Ax.crit = 1.09C7aTff = 0.29 



T \'/Vl03cm-3\'/' 



10 K 



pc, 

(6) 

where Ca is the adiabatic speed of sound in the gas. 
Thus, if MHD waves (or turbulence) contribute even a 
small fraction of the support against gravity (e.g., a few 
percent), ambipolar diffusion removes that support very 
quickly and gravity remains locally unbalanced. Am- 
bipolar diffusion, then, allows thermally supercritical but 
magnetically subcritical masses to begin to separate out 
as fragments (or cores) . Eventually, the mass-to-flux ra- 
tio of the central flux tube of a cloud exceeds the crit- 
ical value for collapse and dynamical contraction (but 
not free fall) ensues. This was the conceptual origin of 
the new theory of fragmentation and star formation sug- 
gested by Mouschovias (1987a; see also 1991a) and put 
on a rigorous footing in subsequent years (e.g., see review 
by Mouschovias 1996 and references therein). 

Statements frequently encountered in the literature, 
to the effect that the ambipolar-diffusion theory ignores 
MHD waves (or turbulence), are thus not based on the 
record. It is precisely because of the partial support 
that such disturbances provide against gravity and their 
consequent decay by ambipolar diffusion that fragmen- 
tation (or core formation) is initiated in this theory. It 
is noteworthy that, the greater the initial contribution 
of such disturbances to the support against gravity, the 
more rapid the ambipolar-diffusion-initiated fragmenta- 
tion process becomes. 

An important prediction of the detailed numerical sim- 
ulations was that the central mass-to-flux ratio, after 
it exceeds the critical value given by equation ([3|) and 
dynamical contraction (but not free fall) sets in, does 
not continue to increase indefinitely. It asymptotes to a 
value typically 2-3 times greater than the critical one 
(e.g., Fiedler & Mouschovias 1993, Fig. 9b; Ciolek & 
Mouschovias 1994, Figs. 2e and 4e). The physical rea- 
son for this is that the contraction is rapid (the infall 
acceleration is typically 30% that of gravity) , so that the 
magnetic flux is essentially trapped inside the contracting 
fragment. Once such gravitational contraction begins, 
unlike the early, ambipolar-diffusion-controUed phase of 
core formation in which the magnitude of the magnetic 
field does not typically increase by more than 30%, the 
magnetic field increases with gas density as B oa p'^ with 
K « 1/2. 

Critics of the ambipolar-diffusion theory often refer to 
the large set of calculations that developed and refined 
the theory as "static". Far from that being the case, 
these calculations have demonstrated that, even if one 
begins with a model cloud that is in an initial exact 



equilibrium state and would therefore remain in such a 
state indefinitely if it were not for the presence of am- 
bipolar diffusion, the fragmentation is relatively rapid 
(more precisely, it occurs on a timescale equal to 1/2 the 
value given by eq. [1]) and, moreover, these calculations 
are the only ones that have followed the subsequent dy- 
namical phase of protostar formation, up to densities of 
about 10^^ cm~3, accounting for the zoo of phenomena 
relevant to the seven-fluid physical system (neutral par- 
ticles, atomic and molecular ions, electrons, negatively- 
charged, positively-charged and neutral grains) - e.g., 
see Tassis and Mouschovias (2007a, b, c). Mindful of the 
possibility that molecular clouds are not quiescent ob- 
jects and are not necessarily found in quiescent environ- 
ments, Mouschovias (1987a, fig. on p. 486; or 1991a, Fig. 
1) considered both subAlfvenic as well as super Alfvenic 
contraction. Yet in both cases, ambipolar diffusion (and 
magnetic braking) plays a crucial role in protostar for- 
mation. 

It was recognized early on in the development of the 
ambipolar-diffusion theory of fragmentation that, once 
a flux tube in a molecular cloud becomes magnetically 
supercritical, it may break up along its length into 
thermally supercritical but magnetically sw^critical frag- 
ments (Mouschovias 1991c, §2.4). In fact, the factor 
by which the density of such fragments must increase 
in order for balance of forces along field lines to be re- 
established and a critical mass-to-flux ratio to be re- 
acquired was calculated analytically. Moreover, the effect 
that such fragmentation has on the B — p relation was 
also determined. It is therefore of utmost importance to 
keep in mind these effects if one desires to make a state- 
ment about the mass-to-flux ratio of a single fragment 
relative to that of the parent cloud. In other words, if a 
magnetically critical or supercritical flux tube breaks up 
along its length, the entire mass in the flux tube must be 
counted if one desires to obtain the relevant mass-to-flux 
ratio to compare with the theoretical prediction. 

The theory has three significant dimensionless free pa- 
rameters, whose typical values are obtained from obser- 
vations (Fiedler & Mouschovias 1992). Accounting for 
rotation introduces a fourth free parameter, which is es- 
sentially a measure of the moment of inertia of the cloud 
(or fragment) relative to that of a region in the surround- 
ing medium (or envelope) of volume comparable to that 
of the cloud (or fragment) (Basu & Mouschovias 1994). 
A simulation is carried out using those typical values of 
the free parameters as input. The calculation predicts 
the spatial dependence and time evolution of the phys- 
ical quantities of the forming protostellar fragment and 
of the parent cloud (e.g., profiles of the density, mag- 
netic field, angular velocity, infall velocity as functions of 
time). In the calculations without rotation, the three free 
parameters are: (1) the initial central mass-to-flux ratio 
/ico in units of its critical value for collape; (2) the ratio 
of the initial free-fall time and the neutral-ion collision 
time, Vfifi] and (3) the exponent k in the parametriza- 
tion of the ion density tij in terms of the neutral density 
n^, k — dhmi/dlnn^. If detailed grain chemistry is used 
to determine the degree of ionization at each stage of 
the evolution, the parameter k is replaced by a quantity 
which is essentially the radius of a typical grain particle 
(Ciolek & Mouschovias 1993). Once the results from the 
fiducial run are at hand, a parameter study is carried 



out. It consists of at least two additional runs for each 
of the free parameters, using values at the endpoints of 
the range suggested or allowed by observational uncer- 
tainties. Since the mass-to-flux ratio of the parent cloud 
(or envelope) is a dimensionless free parameter in the 
theory, the theory neither requires nor predicts that this 
quantity have any specific value, e.g., be subcritical. The 
fact that, in any cloud in which a positive detection of 
the magnetic field component along the line of sight using 
Zeeman observations has been obtained, implies that the 
cloud's mass-to-flux ratio (given the geometrical uncer- 
tainties in estimating it) is consistent with a subcritical 
value is just that, observational input to the theory, not 
a requirement or a prediction of the theory. 

Early observational work using the Zeeman effect in the 
21-cm line of HI, focused on testing the relation between 
the magnetic field strength and the gas density predicted 
by spherical, isotropic contraction, i.e., B oc p^'^ (e.g., 
Verschuur 1970, 1971). It was concluded that the obser- 
vations were consistent with that relation. However, a 
more careful examination of the data showed that that 
was not the case (Mouschovias 1978, Fig. 1). Instead, 
the data was consistent with the prediction of the early, 
self-consistent calculations of the equilibrium and con- 
traction of self-gravitating, isothermal, magnetic clouds 
embedded in a hot and tenuous external medium. Sub- 
sequent observations of the Zeeman effect in OH (which 
samples greater densities than HI) and in OH and H2O 
masers (which sample even greater densities and reveal 
physical conditions in regions of active star formation) 
were found to be as predicted by the magnetic calcula- 
tions (Fiebig & Guesten 1989; Mouschovias 1996, Fig. 
lb; Crutcher 1999). 

More recently, observational attention has shifted from 
the B — p relation to the mass-to-flux ratio, which, for 
a cloud as a whole is invaluable input to the ambipolar 
diffusion theory and, for fragments (or cores) can poten- 
tially be an important test of the theoretical prediction 
that ambipolar diffusion leads to an increase of the mass- 
to-flux ratio of a cloud's central flux tubes relative to the 
mass-to-flux ratio of the envelope. 

2.2. Recent OH Zeeman Observations 

CHT have recently combined existing detections of the 
OH Zeeman effect in four molecular cloud cores with new 
observations of the same effect in the regions surrounding 
these cores, in an effort to get a handle on the variation 
of the mass-to-flux ratio from the envelope to the core of 
each cloud. The authors regard this as a definitive test 
of the ambipolar diffusion theory. The four clouds are 
L1448, B2 17-2, L 1544, and Bl. In the region surround- 
ing each core in each of the four clouds, new observations 
were undertaken at four locations about each core (i.e., 
in the envelope of each core) . At some locations in each 
envelope, positive detections were made while at oth- 
ers no detection could be achieved. For each envelope, 
the authors average the algebraic values of the four mea- 
surements of the line-of-sight magnetic field, regardless 
of whether they were detections or nondetections, and 
assign the resulting value to the average magnetic field 
strength of the envelope - the assumption made here 
is that this value appropriately characterizes the mean 
magnetic field of the envelope. Using this field, CHT 
proceed to calculate what they regard as the magnetic 



flux of the envelope, which, combined with the flux in 
the core, is used to obtain the quantity R defined by 

j-f -'corc^ *^corc/ -^corc /„\ 

^^~T AT? Tr ■ ^'> 

(The quantity / is the peak intensity of the spectral line 
in degrees K, AV^ is the FWHM in km s^^, and B is 
the line-of-sight magnetic field strength in microgauss.) 
They then average again the values of R over the four 
clouds, and obtain an overall average value of R. The 
claim in the paper is that the derived value of R is Se- 
away from 1 and that this contradicts a prediction of 
the ambipolar diffusion theory and, therefore, invalidates 
the theory, since the latter (according to CHT) predicts 
R>1. 

We show quantitatively below that this claim, as stated 
by CHT, is not supported by their data; it is an artifact of 
the analysis and interpretation of the data, which suffers 
in three major ways: 

1. The treatment of the propagation of observational 
uncertainties systematically underestimates the uncer- 
tainties of the combined result at every step. 

2. Nondetections of magnetic fields are treated as if 
they were detections, and upper limits are not quoted as 
appropriate - in fact, based on the raw data presented 
in CHT, measurements of R are not achievable in any of 
the clouds studied; rather, upper limits should have been 
calculated and quoted. 

3. The comparison with ambipolar-diffusion calcula- 
tions is itself flawed in three important respects: 

(a) The theoretical simulations chosen for this compar- 
ison use as input a very different geometry of the field 
lines threading the cloud from the geometry of the field 
lines that the observations presumably reveal. 

(b) The magnetic flux threading each cloud is incor- 
rectly calculated by taking the arithmetic average of the 
four algebraic values of B in each envelope. 

(c) The possible breakup of the mass in the central flux 
tube of each cloud into thermally supercritical but mag- 
netically subcritical fragments (or cores), an important 
possible effect in the theory (Mouschovias 1991c), is men- 
tioned by the authors but not taken into consideration 
in their analysis and/or arguments. 

If these effects are considered, it follows that a variety 
of i?— values are possible. We discuss each of these points 
below. 

3. COMBINING UNCERTAINTIES AND PROPAGATING 
ERRORS 

CHT quote values for the results of observations of the 
magnetic field at four different positions in the envelope 
of each cloud, together with associated uncertainties, de- 
rived as in Troland & Crutcher (2008). In most cases, 
the quoted values correspond to nondetections. We dis- 
cuss this issue in detail in the next section; here, we refer 
to detections and nondetections alike as "observations" . 

For each observation of the envelope's line-of-sight 
magnetic field Bj , CHT quote an associated uncertainty 
<Tj, which we assume to be Gaussian. These values are 
given for each cloud in Table [1] The authors then cal- 
culate a mean value of the line-of-sight magnetic field in 
the envelope of each cloud. 
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TABLE 1 

Magnetic Fields and Errors (in microgauss) in Four 
Cloud Envelopes (data from CHT). 



Cloud 


Bl±(Tl 


B2 ± (72 


B3±CT3 


B4±cr4 


L1448CO 


-9 ±13 


-11 ±6 


-7 ±7 


14 ±8 


B217-2 


-13 ±9 


5±6 


6±8 


9 ±13 


L1544 


-3 ±4 


-1±4 


22 ±6 


2 ±10 


Bl 


-16 ±6 


0±7 


-3 ±6 


-10 ±5 



and an uncertainty on the mean, (Tnioan, using Gauss' 
formula of error propagation on equation ([8]) , which gives 
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The results derived using this formula are given in the 
second column of Table [51 they are the same as those 
quoted in CHT (see their Table 1, col. 4). This analysis 
is a standard treatment of errors when combining several 
measurements of the same quantity. If CHT had mesured 
the magnetic field at the same envelope location in four 
different instances, their measurements could be com- 
bined using this procedure (although, even in this case, 
a weighted average and its associated error propagation 
formula would be more appropriate since the measure- 
ment uncertainties are not equal for every observation). 
Similarly, if there were an a priori reason to be certain 
that the envelope magnetic field is uniform and the four 
observations in each cloud are sampling the same value 
of the field, such treatment would be appropriate. How- 
ever, this is not the case here, and the treatment of errors 
in CHT yields misleading results. We discuss the reasons 
below. 

First, it cannot be overemphasized that Gauss' for- 
mula for error propagation, and hence equation ([9|), is 
not universally applicable. There are two reasons for 
which applicability might fail. The first stems from the 
fact that Gauss' formula is based on a Taylor expansion 
(e.g., see. Wall & Jenkins 2003). As such, it is valid only 
if the errors being propagated are small, or if the equa- 
tion through which the errors are propagated is linear 
(and thus no second-order terms enter the calculation). 
This condition does indeed apply in this case: although 
the errors are by no means small, the average is a lin- 
ear operation and no second-order terms enter the error 
propagation formula. The second reason for which the 
applicability of Gauss' formula might fail is if the ob- 
servations being combined are not measuring the same 
intrinsic quantity, but rather a distribution of values. 
Then, one's ability to constrain the mean, especially in 
the small-number-statistics regime, is not simply depen- 
dent on one's ability to take individual measurements; 
the possibility of having finite distribution spreads needs 
to be explicitly accounted for. 

The magnetic field in the cloud envelope is not known 
a priori to have a unique uniform value everywhere in 
the envelope. In fact, the data seem to suggest exactly 
the opposite (e.g., compare observations 1 and 2 in cloud 
L1544 with observation 3 in the same cloud - see Table 
[T|) . Equation ^ is therefore not the appropriate estimate 
for the uncertainty on the mean magnetic field for each 



cloud envelope. What should one use instead? 

The spread in the observations is the convolution of 
the observational uncertainty and the intrinsic spread of 
the distribution of magnetic field values across the enve- 
lope. The sample variance can give an estimate of this 
convolved spread; its square root (the sample standard 
deviation, an unbiased estimator of the population stan- 
dard deviation a) is given by 
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It is well known (e.g., see Wall & Jenkins 2003; Lyons 
1992) that the average, given by equation ([H]), is dis- 
tributed as a Gaussian about the true value of the mean, 
with variance u"^ /N where N is the number of observa- 
tions. In our case, A^ = 4 and the standard deviation of 
the average, which is the uncertainty on the mean, is a /2. 
For all clouds, the best-guess estimate of a, which is as, 
is > twice the uncertainty quoted on the mean (note, for 
example, the case of L1544 where the uncertainty quoted 
by CHT on the mean is 3^G while (Ts/2 — 6/iG). This 
is an empirical proof that the uncertainties quoted by 
CHT are underestimated, even if one were to ignore the 
conceptual point made above. 

Even (75/2, however, underestimates the true uncer- 
tainty on the mean. The uncertainty on the mean is (t/2, 
not <Ts/2. The sample standard deviation cr| is only an 
estimate of the population variance u^ , and is distributed 
as a^^ l(N — 1), where x^ is a X"S<luare variable with 
A^— 1 degrees of freedom (e.g., see Wall & Jenkins 2003). 
In practical terms, this means that, in the low-statistics 
regime we are considering, ag most frequently underesti- 
mates cr, and in about 10% of all instances does so by at 
least a factor of 2. In the latter case, the appropriate un- 
certainty on the mean is closer to as rather than as/'i-- 
How do we then take into account not only that there 
is a finite population spread, but also that this popula- 
tion spread is essentially unknown, due to our inability 
to constrain it very well with the very limited number of 
available observations? 

We must use an analysis that allows not only for finite 
spreads, but also for a variety of such spreads, as long 
as they are allowed by the data. Allowing for a range of 
finite spreads consistent with the data increases the range 
of mean values that can explain the data, and in this 
way increases the uncertainty on the mean. Conversely, 
ignoring the possibility that there is true variation of 
mostly unknown extent of the B\os from point to point 
in the envelope artificially decreases the uncertainty in 
the mean i?ios value. 

The most straightforward way to correctly combine 
these observations while allowing for finite spatial spread 
of Bios is through a likehhood analysis (e.g., see Wall & 
Jenkins 2003; Lyons 1992; Lee 2004). In a likehhood 
analysis, the data can be assumed to be derived from 
a parent distribution of values (in our case, we assume 
that values of i?ios in the envelope follow a Gaussian dis- 
tribution with mean equal to Bq and intrinsic spread 
fJo)- This distribution is then "sampled" with A^ mea- 
surements Bj , each carrying a (Gaussian) uncertainty of 
measurement aj. ^^From the data, we then calculate the 
joint probability distribution of the parameters BQ,ao 
which are consistent with observations (the likelihood). 



TABLE 2 
Derived mean envelope magnetic fields and uncertainties. 

Two left COLUMNS: CHT TREATMENT; THIRD COLUMN: SAMPLE 

STANDARD DEVIATION; FOURTH COLUMN: LIKELIHOOD ANALYSIS 

RESULTS (this WORK). 
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In our case, the likelihood can be shown to be (e.g., see 
Venters & Pavlidou 2007, and appendix R|) 
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Any parameters that are not of direct interest to us (such 
as (To in our case), can then be integrated out of the like- 
lihood. In this way, we can derive the probability dis- 
tribution of the parameter of interest (Bq in our case) 
while still allowing for all possible values in uq. The in- 
tegrated likelihood is called the marginalized likelihood^ 
Cm] this probability distribution can then be used to de- 
rive confidence intervals and upper limits where appro- 
priate. A detailed discussion of the likelihood approach 
to error analysis, and the use of the marginalized likeli- 
hood to derive confidence intervals and upper limits, is 
given in Appendix El An example of the (unnormalizcd) 
marginalized likelihood for the case of L1448CO is shown 
in Figure[T] The marginalized likelihood is derived by nu- 
merically integrating equation (|11[) over co for different 
values of Bq, and is shown as a solid line; the location 
of the maximum-likelihood estimate for the mean Bq is 
indicated by the dash-dot line. 

The maximum-likelihood estimates and associated un- 
certainties of Bo are shown in Table [21 It is clear that 
these are systematically greater than the errors produced 
by (the nonapplicable in this case) equation ([9]) , and they 
are generally greater than (Ts/2. The value of £mdefining 
the 1(7 uncertainties in the case of L1448 and the associ- 
ated -Bo- and i?o+ sue marked with the dashed lines in 
Figure [H For comparison we show, with the shaded box, 
the 1(7 spread of the values of Bq that CHT quote for the 
same object, based on the same data. 

Besides the significant underestimate of uncertainties 
in the mean envelope magnetic field quoted in CHT, the 
second striking feature seen in Figure [1] is that the "mea- 
sured" mean envelope magnetic field is very much consis- 
tent with zero (even using the CHT uncertainties, zero 
is marginally more than one a away from their "best- 
guess" mean), at least in the case of L1448CO. It is also 
clear from Table [2] that, the maximum-likelihood values 
of Bq for all clouds are no more than 2a away from zero. 
As a result, the observational data presented in CHT are 
inconsistent with a mean value of the envelope magnetic 



field equivalent to a "detection" in any of the observed 
clouds (note however that in certain individual locations 
within cloud envelopes detections have been achieved; 
the significance of these results is nevertheless diluted by 
the rest of nondetections, and the mean is never more 
than 2ct away from zero). Given this fact, one would ex- 
pect to see appropriate upper limits quoted for the mean 
magnetic field in the envelope for each of these clouds. 
Such limits are not given; instead, the "best-guess" val- 
ues and their uncertainties are treated in exactly the 
same way that one should treat only detections. This 
makes the interpretation of the data presented in CHT 
problematic. We discuss this issue in the following sec- 
tion. 

4. UPPER LIMITS VERSUS DETECTIONS 

A "detection" of a magnetic field is an observation with 
a high enough signal-to-noise ratio that the measured 
magnetic field is statistically significantly different from 
zero. In other words, when a detection occurs, we are 
confident at some level (99.7% for a nominal 3(7 detec- 
tion) that the magnetic field measurement is not spuri- 
ous (a result of noise). If a detection has indeed occured, 
then the magnitude and sign of the magnetic field can 
be quoted, together with associated error bars (indicat- 
ing the Ict uncertainty in the magnetic field value). If a 
detection has not occured, a value of the magnetic field 
is meaningless and is therefore not quoted; what is rou- 
tinely quoted instead is an upper limit (a value of the 
magnetic field magnitude which, if exceeded in the ob- 
served system, would lead to a detection). 

In fields of science such as particle physics detections 
and upper limits are routinely quoted at the 5a level; 
even in the literature of Zeeman observations of inter- 
stellar magnetic fields it is common practice to adopt 
as detections 3cr measurements, and accordingly quote 
3(7 upper limits where appropriate (e.g. Crutcher et al. 
1993). However, for consistency with CHT and what 
they consider as a tentative detection, we adopt 2a limits 
here. Formally, then, the 2a upper limit on the magnetic 
field is the value of the magnetic field for which we are 
confident at the 95.4% level that it is not exceeded in the 
observed system. 

Even if we use the (significatly underestimated, as we 
discussed in the previous section) overall uncertainties of 
CHT quoted in Tabled and treat them as perfectly gaus- 
sian errors (which, as we have seen, is not the case), 3 of 
the 4 clouds do not have a statistically significant (> 2a) 
measurement of the mean envelope magnetic field, while 
the 4th cloud (Bl) has a marginal, 2a measurement of 
the mean; even the Bl mean envelope magnetic field be- 
comes consistent with zero if the correct uncertainties are 
used. As CHT point out, this is not necessarily a prob- 
lem: useful limits on the ratio R can be derived using up- 
per limits for the envelope magnetic field, provided that 
measurements of the core magnetic field exist. However 
such limits, either for the envelope magnetic field or for 
R, are not calculated by the authors. Instead, nondetec- 
tions are treated as detections, and the authors repeat- 
edly refer to their "measurements" of R, R' , and their 
mean. Such measurements are not, in fact, possible with 
the data presented in their paper. The ratio R has not 
been measured. However, upper limits can be placed on 
its value although, as we explain below, these limits are 
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TABLE 3 

2(T Upper Limits on the Envelope Magnetic Field (in fiG) 

AND on R. 
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Fig. 1. — Marginalized likelihood for the case of L1448CO (not 
normalized, heavy black line). The location of the maximum- 
likelihood value of Bo is shown with the dot-dashed line. The 
dashed lines indicate the la equilikelihood level, and the locations 
of the 1(T limits of Bq. The shaded box shows the allowed Icr values 
for the mean envelope magnetic field for the same cloud according 
to CHT. Dotted lines indicate the 2cr upper limit on the absolute 
value of the mean envelope magnetic field. 



not very strong for any of the observed clouds. 

We note that the proper treatment of nondetections 
can in principle severely constrain theories of interstellar 
medium evolution and core formation, if the upper limits 
placed are strong enough (it could, for example, be the 
case that the measured 3cr upper limit for R is 10~^! - 
that would certainly be implying a very weak envelope 
magnetic field and a very puzzling situation for any cur- 
rent core formation theory). However, the magnetic field 
upper limits have to be correctly calculated and treated 
as such - and they do have to be strong enough, which 
is not the case at all for the data presented in CHT. 

Upper limits for the (absolute value of the) mean mag- 
nitude of the line-of-sight magnetic field combining all 
measurements of the envelope magnetic field can be de- 
rived starting from the marginalized likelihood of the pre- 
vious section. The procedure is discussed in detail in 
Appendix |^ Utilizing the marginalized likelihood pro- 
vides a very straightforward way to account for and cor- 
rectly treat the asymmetric nature of the error bars in the 
combined observational results for the envelope magnetic 
field, the long non-gaussian tails, and the exact nature 
of the question we are trying to answer. 

The 2(7 upper limits for the envelope magnetic field 
calculated this way are given in Table [3] The values of 
I Bo I and — |i?o| which include between them a fractional 
area of 2a of the marginalized likelihood for L1448CO 
are marked with dotted lines in Figure [T] Note that the 
marginalized likelihood at — |-Bo| and at |i3o| need not 
have (and does not have in this case) the same value; 
this is frequently encountered when the robust definition 
of an absolute-value mean envelope B-field upper limit 
is used. 

Since no measurements of the mean envelope magnetic 
field |i?cnv| have been achieved, no meaningful values of 



Cloud 


|Senv|(<2(T) \R 


(<2a) 






(/.G) 






L1448CO 




27 




2.0 


B217-2 




22 




2.9 


LI 544 




29 




5.0 


Bl 




20 




1.1 



the ratio R can be inferred. We can, however, use the 
derived upper limits to similarly place upper limits on 
the values of R. To treat the problem properly, and cor- 
rectly propagate uncertainties on the measured or lim- 
ited observables to the derived quantity R, we cannot use 
equation ([9|), since the individual uncertainties involved 
are large while at the same time the equation which de- 
fines R (equation [7]) is nonlinear. Instead, we do a full 
Monte- Carlo calculation to properly derive the probabil- 
ity distribution for the values of R as follows. We repeat 
the following experiment 10^ times: we draw /core, ^cnv, 
AT4orc, Acoro and Score from gaussian distributions with 
mean and spread equal to the measurement and uncer- 
tainty quoted in CHT (this is equivalent to assuming all 
errors in these quantities to be Gaussian); we draw a 
mean value of i3onv from the marginalized likelihood of 
the previous section; we combine all the "mock observa- 
tions" of these numbers to produce one value of R. We 
use the 10^ values of R produced in this way to numer- 
ically calculate the probability distribution for R. We 
then calculate the 2a upper limit on \R\ by requiring 
that the fractional integral of this distribution between 
—R and R be 95.4%. The 2a upper limits for \R\ are 
given in Table [31 Unfortunately, none of these limits are 
very strong. Although this is an interesting thought ex- 
periment, much more and better quality data are required 
before any statement can be made with regard to compar- 
ison with any theory. 

In the data presented by CHT, there is a single statisti- 
cally significant, at the 3a level, detection of the envelope 
magnetic field - one of the four observations of the L1544 
envelope magnetic field yields a value of 22 ± 6 /iC This 
measurement is consistent with the 29 fiG 2a upper limit 
that one derives if all four observations of the envelope 
magnetic field are combined for the particular cloud. If 
we use this value with its associated uncertainty to calcu- 
late an R for the case of L1544, we obtain 3.62 ± 1.29 - a 
marginal detection of the ratio R, with a value of a few. 
Again, this value is consistent with the derived overall 
3(7 upper limit of 5.0 obtained for L1544 by combining 
all observations of the magnetic field. 

5. COMPARISON WITH THE AMBIPOLAR-DIFFUSION 
THEORY 

The main conclusion of the CHT paper is that their 
measured values of R are significantly smaller than 
1, which is in contradiction with a prediction of the 
ambipolar-diffusion theory. We have shown above that 
there are, in fact, no measurements of R possible with 
the CHT data; only weak upper limits can be placed, 
which have values of a few. This, however, is only part 
of the problem with their strongly stated conclusions. 




Fig. 2. — Schematic diagram of a star (e.g., Sun) that has a 
dipolar magnetic field and is observed along lines of sight parallel 
to its equatorial plane. To calculate the magnetic flux threading 
the star, only the mean value of the magnetic field in either the 
northern or in the southern hemisphere is needed. If both values 
are averaged algebraically, the erroneous fiux value of zero will be 
obtained. 



The second error steins from the manner in which they 
calculate the magnetic flux of each cloud, even if their en- 
velope magnetic fields and associated errors were to be 
exactly as quoted in their paper. They add algebraically 
the four values of the "measured" line-of-sight magnetic 
field strength and divide by four to obtain the mean 
field of the envelope. (Then they multiply this value 
by the area of each envelope in the plane of the sky to 
get the magnetic flux.) This procedure badly underesti- 
mates the envelope's magnetic flux in at least three of the 
four cases, in which the envelope "measurements" show 
fields with opposite algebraic signs at different positions. 
These reversals in the direction of the magnetic field of 
the envelope, if real, would imply a bent magnetic flux 
tube threading the cloud. In such a case, only one alge- 
braic sign of the magnetic field (the one corresponding to 
the greatest absolute values) should be considered in es- 
timating the magnetic flux of the envelope. An example 
will clarify the point. Suppose that the dominant com- 
ponent of the magnetic field of a star, such as the Sun, 
is dipolar in nature and suppose that one wants to esti- 
mate the magnetic fiux threading the star by measuring 
the magnetic field on the star's surface. We assume, for 
simplicity, that the line-of-sight from the observer to the 
star lies in the star's equatorial plane, as shown in Figure 
[21 If this observer did what CHT do, he would conclude 
that the magnetic fiux of the star is exactly zero. That, 
of course, is far from being the case. The magnetic fiux 
threading the star is obtained from the mean value of 
the magnetic field of either the sourthcrn or the north- 
ern hemisphere of the star (but not from the algebraic 
sum of the two) by multiplication with the area of the 
equatorial plane of the star. 

Moreover, although CHT believe that their data re- 
veal reversals in the direction of the magnetic field within 
the clouds they observe, they nevertheless compare their 
results with calculations that follow the formation and 
evolution of fragments in molecular clouds under the as- 
sumption that the parent clouds' magnetic field lines 
are initially essentially straight and parallel. This ide- 
alization in the theoretical calculations renders a mathe- 
matically complicated multifiuid, nonidcal MHD system 
tractable while capturing all the essential physics of the 
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Fig. 3. — Schematic diagram (from Mouschovias and Morton 
1985b) of a deformed fiux tube that has fragmented along its length 
in a molecular cloud. The deformation can be caused by the rel- 
ative motion of the fragments. The fragments need not be mag- 
netically connected for their relative motion to cause significant 
deformations of the field lines. The dashed lines represent different 
lines of sight, whose significance is explained in the text. 




Fig. 4. — Schematic diagram of the deformation of the field lines 
threading a cloud caused by its motion relative to the surrounding 
medium. The cloud is shown, for simplicity, to contain only one 
fragment (or core), in the neighborhood of which the hourglass 
shape of the field lines had been established during core formation 
but affected by the cloud's motion. The dashed lines represent 
different lines of sight, whose significance is explained in the text. 



core formation and evolution problem. However, nobody 
has ever suggested that in a real cloud the magnetic field 
lines will be essentially straight and parallel, except for 
the local hour-glass distortions associated with the com- 
pression of the field lines during gravitational core for- 
mation. In fact, Mouschovias & Morton (1985b, Fig. 
13) had sketched what they regarded as a more realistic 
field geometry in a molecular cloud in which there are 
several (in that case four) magnetically connected frag- 
ments. That figure is reproduced here as Figure [3l This 



configuration can result from relative motion of the frag- 
ments (labeled A, B, C, and D in the figure) within the 
cloud. The motion of a cloud as a whole relative to the 
intercloud medium will also bend the magnetic field lines 
in an almost U shape, as sketched in Figure 4. One can 
easily visualize lines of sight in Figures [3] and |4] along 
which, if one were to observe -Bios 7 one would detect re- 
versals and/or essentially vanishing -Bios- In other words, 
observations that can potentially reveal the geometry of 
the field lines can and should be used as input to make 
refinements in the theory and/or build a particular model 
for the cloud observed (as, for example, done in the case 
of Bl by Crutcher et al. 1994, and in the case of L1544 
by Ciolek & Basu 2000); such observations by no means 
invalidate the theory, which afterall has produced dozens 
of predictions many of which have been confirmed by ob- 
servations and none of which has been contradicted by 
any observation. 

Figures [3] and [5] also demonstrate that the assumption 
by CHT, namely, that cos9, the geometrical correction 
accounting for the fact that only the line-of-sight mag- 
netic field is measurable by Zeeman observations, is the 
same in the core and everywhere in the envelope, is not 
in general valid. 

The lines-of-sight labeled AA', BB', CC, DD' and EE' 
in Figure 4 are chosen so as to reveal qualitatively sim- 
ilar behavior of the core and envelope -Bios as the corre- 
spondingly labeled lines-of-sight in Figure [31 In relation 
to Figure m if an observer's line-of-sight is the line DD', 
he'll detect essentially the full strength of the magnetic 
field in the nearest side of the cloud's envelope. Along 
CC, he will detect a weaker -Bios in the core, although the 
full strength of the core's magnetic field is clearly greater 
than that of the envelope, as revealed by the compressed 
field lines. By contrast, along the EE' line-of-sight no 
envelope -Bios will be detected; and along AA', almost 
the full strength of the core's magnetic field will be mea- 
sured, but only a fraction of the envelope's magnetic field 
strength will be detected. It is clear from this illustrutive 
example, which actually may not be very different from 
the field line geometry in a real cloud (e.g., see Brogan 
et al. 1999), that unless one is mindful of geometric ef- 
fects, incorrect conclusions will be reached on the true 
magnitude of the magnetic field and, therefore, of the 
mass-to- fiux ratio in cloud cores and envelopes. When 
one recalls that the relevant column density for calcu- 
lating the mass-to-flux ratio is the one along a mangetic 
flux tube, not the one measured along an arbitrary line- 
of-sight, which may be at a very large angle with respect 
to the true magnetic field direction, one realizes that geo- 
metric effects cannot be ignored in attempts to measure 
or estimate the mass-to-flux ratio of a cloud or a core 
(e.g., see Shu et al. 1999). Simply taking the ratio of 
the two does not eliminate the need to account for the 
geometric effects. 

6. SUMMARY AND CONCLUSION 

We have examined critically the analysis and conclu- 
sions of recent OH Zeeman observations of four molecular 
cloud envelopes and cores by Crutcher et al. (2008), and 
we have shown that: 



1. The error analysis of measurements of the envelope 
magnetic field is flawed, and the uncertainties on 
the observed mean value of the field are systemat- 
ically underestimated in every cloud. 

2. Although the mean value of the envelope magnetic 
field is not statistically significantly different from 
zero in any cloud, the authors treat these values 
as if they were detections, and they do not give 
appropriate upper limits. 

3. If the error analysis is performed correctly, allowing 
possible spatial variations of the line-of-sight mag- 
netic field (suggested by the data), and upper hmits 
are appropriately calculated for the mean magnetic 
field in the envelopes, these upper limits are very 
weak: the mean magnetic field in the envelopes is 
constrained to be smaller than 20 — 30 jjlG at the 
2(7 significance level. 

4. If corresponding 2a upper limits are calculated for 
the value of the ratio -R, which is given by equation 
(7), they are also weak: R is constrained to be 
smaller than a few. For no cloud is the upper limit 
smaller than 1. -R is never statistically significantly 
different from zero, and has not been "measured" 
for any cloud. 

5. The magnetic fiuxes of the clouds were calculated 
incorrectly, by averaging algebraically the (pre- 
sumed) values of the line-of-sight magnetic field 
-Bios over a cloud's area projected on the plane of 
the sky, despite claimed reversals in the field direc- 
tion not along one and the same line-of-sight, but 
from one line-of-sight to another. 

6. The possible gravitational breakup (of the mass 
along the flux tubes threading the cores) into 
thermally supercritical but magnetically subcriti- 
cal fragments, an effect of the ambipolar-diffusion 
theory studied early on (Mouschovias 1991c), was 
ignored by the authors. Although it is not clear 
whether this effect has occurred in the four clouds 
studied, it is nevertheless important to consider it 
properly before bold conclusions are drawn on the 
validitiy of the ambipolar- diffusion theory. 

7. The ambipolar-diffusion theory can explain values 
of -R smaller than 1, especially if three-dimensional 
geometrical effects are properly accounted for. 
Although such effects are indeed missing from 
the published (axisymmetric) calculations against 
which CHT compare their observations, this is 
hardly new information. Focusing on effects that 
are known a priori to be missing from specific the- 
oretical calculations is not a proper way to test 
the underlying theory or to make observational 
progress toward distinguishing between alternative 
theories. Such observations, when reliable, would 
provide useful input to the ambipolar-diffusion the- 
ory. 
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APPENDIX 
LIKELIHOOD APPROACH TO ERROR ANALYSIS 

In the classical sense, the likelihood represents the probability to have observed the data given some parent probability 
distribution. In the Bayesian sense, the likelihood is also proportional to the probability of the parent distribution 
to be described by a specific set of parameters (mean, variance) given the data. In a likelihood analysis such as the 
one presented here, the likelihood (properly normalized) is the probability of the set of parameters we are trying to 
estimate (mean, variance), as all Bayesian priors are assumed to be flat ^. 

For simplicity, we assume that values of Bios in the envelope follow a Gaussian distribution with mean equal 
to Bq and intrinsic spread cto ^- If we could observe magnetic fields with zero observational error, the proba- 
bility to observe a value B at some specific envelope location (the likelihood I of a single observation) would be 
I — exp [—{B — _Bo)^/2(To] /V^o'o- Since we have observational error, complications arise. At any specific envelope 
location, there is a probability exp [—{B — Bo)^/2ctq] /-\/27r(To for the magnetic field to have a true value Bq- If the 
(Gaussian) observational uncertainty at this same location is CTj, then the probability to observe a value Bj of the field 
given that its true value is B is the latter probabihty multiphed by a factor exp [—{B — i3j)^/2cr|] /^/2TTaj. But this 
is not the only way we could get an observed field value Bj. The true value of the field might be B' , and we would 
still have a finite probability to observe Bj (same reasoning as above, with B substituted by B'). To get the total 
probability for a single observation of Bj , we need to integrate over all possible "true" values of the magnetic field at 
a single location. Doing this, the likelihood for a single observation Bj with observational uncertainty CTj becomes 

,,R , /--^ exp [-{B - B,)y{2a])] exp [-{B - Bo)y{2a'„)] 

lj[Bo,ao) = / dB -= y= . (Al) 

J-oo y Ziraj y Zttctq 

The likelihood L for N observations of Bj with individual gaussian uncertainties aj to come from a Gaussian intrinsic 
probability distribution with mean Bq and spread ctq is simply the product of the individual likelihoods, C = Y[j=i h 
which, after performing the integration in equation (jAip and some algebraic manipulations, yields (e.g., see Venters 
& Pavlidou 2007) 



N 



C{Bo,ao) = Y[ I ^^^P 



, j=i \i (Jn -r a 



1 Y^ {Bj - Bq) 



2 ■^ <Jn + (J 



2 



(A2) 



Note that for ctq = 0, the value of Bq that maximizes equation (|lip is exactly the weighted average, as it should be: 
the weighted average formula is the maximum-likelihood estimator of the mean when the spread comes purely from 
observational uncertainty and the observational uncertainties are unequal. 

The likelihood of equation ([TT|) is a powerful tool that allows us to calculate the joint probability distribution of 
the mean Bq and the unknown intrinsic population spread ctq. However, if one wants to calculate the magnetic flux 
passing through a surface on which the magnetic field has a probability distribution of values, one finds that the flux 
integral is exactly equal to the product of the distribution mean Bq times the surface area (in this case, modulated by 
cos 6* since only the line-of-sight component of the magnetic field is measured). Although we need to allow for different 
values of ctq to correctly estimate the uncertainty on Bq, the actual value of ctq is not of interest for the purpose of 
this analysis. We can therefore integrate it out of the likelihood (compute the marginalized likelihood of the mean, 
'Cm = / CddQ , which physically corresponds to calculating the probability distribution of the mean for any value of uq 
allowed by the data). The marginalized likelihood gives the probability distribution of Bq allowing for any possible 
(To that can explain the data. The marginalized likelihood is asymmetric, and so will be the error bars. In addition, 
the location of the maximum-likelihood value of Bq docs not have to agree, and generally does not agree, with the 
estimate of the more restrictive equation (jS)) . 

To calculate uncertainties of the mean, we start from the peak of the marginalized likelihood and progressively lower 
the likelihood level, calculating the fractional integral under the curve and between the two values i?o_ and _Bo+ left 
and right of the maximum-likelihood value, which are such that £^(^0-) = C,m{BQ+). We do so until we find the 
Cm value for which this fractional integral is equal to 68.26%. The -I- and — la uncertainties are then defined as 
Bq+ - B^^^c and Bmax£ - -Bq-, respectively. 

Upper limits can also be calculated starting with the marginalied likelihood, as follows. The 2cr upper limit is that 
value of \Bq\ for which the area under the marginalized likelihood curve between — |i?o| and |i?o| is 95.4%. In other 
words, a fractional area equivalent to the gaussian 2a is included between — |-Bo| and |i?o|- Note that, because of the 
different nature of the question (in this case we are interested in fixed values of Bq rather than equal £,„ likelihood 
levels), this is not simply the maximum- likelihood value of Bq plus two times the error. In addition, the marginalized 

^ or, equivalently, flat everywhere where the likelihood has appreciable power; setting upper limits based, for example, on prior obser- 
vational knowledge that the magnetic field in the interstellar medium is less than tens of thousands of Gauss can yield a flat yet proper 
Bayesian prior, if such a property is desired. Here, the results are completely independent of such a choice. 

^ Since we still have to assume a functional form for the probability distribution of -Bios within the envelope, the uncertainty we derive 
for Bq is still underestimated since we are not accounting for the possibility that -Bios is distributed according to some other functional 
form. 
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likelihood has much longer tails than a gaussian, and hence la, 2ct, ... values do not scale linearly (if we desire to keep 
the terms representing the same fractional integrals as in the case of a gaussian). 
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